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In baseband combining , the key element is the time alignment of the baseband 
signals . rhis article describes a closed-loop time-alignment system that estimates 
and adjusts the relative delay between two baseband signals received from two differ- 
ent antennas for the signals to be coherently combined. This system automatically 
determines which signal is advanced and delays it accordingly with a resolution of a 
sample period . The performance of the loop is analyzed , and the analysis is verified 
through simulation. The variance of the delay estimates and the signal-to-noise 
ratio degradation in the simulations agree with the theoretical calculations. 


I. Introduction 

Antenna array combining techniques have been used 
in the Deep Space Network to improve the signal-to-noise 
ratio (SNR) [1]. These techniques include full spectrum, 
baseband, and symbol stream combining. Each technique, 
however, has been applied only to a specific mission, and a 
rigorous comparison between the techniques is missing. To 
have a better understanding of how these techniques per- 
form under the same or different conditions, all techniques 
are being studied and simulated in software. 


This article concentrates on baseband combining where 
signals from two antennas have each had their carrier re- 
moved before going into the combiner. A diagram of base- 
band (BB) combining is shown in Fig. 1. The signals re- 
ceived from the zth antenna have the following form [1]: 


r -(<) = \fiPc, sin (w c ,(< - r.) + 0 C> ) + s/2P D% d(t - r,) 
x sgn [sin (u, c (t — 7 \) -f- 0 JC )] 

x cos(w Ci (i — tv) + e Ct ) + m(t) (1) 

where P c , = P, cos 2 A and P D> = P; sin 2 A, with P { being 
the total signal power of the received signal from the ith 
antenna, and A is the modulation index. The parameters 
u> C| and 6 Ci denote the carrier angular frequency and phase 
of the received signal from the ith antenna, u; 5C and 0 SC 
denote the subcarrier angular frequency and phase, and 
T{ is the relative delay between the reference signal and 
the ith signal. These received signals are analog- to-digital 
(A/D) converted and downconverted to an intermediate 
frequency (IF) and then sampled before going through the 


92 




CPLRS + PAD 



carrier tracking in the model. After carrier removal, ig- 
noring the filtering effects on the data and the squarewave 
subcarrier, the signals have the form 

r' t ( k ) = y Po t cos 4> H d{kT, - r,) 

x sgn [sin ( u sc {kT, - r,) + 9 SC )] + nJ-(A) 

= Si(k) + n[(fc) (2) 

where T s denotes the sampling period, and (f) Ct is the car- 
rier phase error in the Mi carrier-tracking loop. The sig- 
nals r-(£), i — 1,2,..., are the to-be-combined baseband 
signals with different time-varying time delays. To coher- 
ently combine these signals, it is necessary to align them 
in time first. 

To align two signals in time, first find the relative de- 
lay between the two signals. A standard procedure for 
finding a relative delay between two signals is to correlate 
them [2] since the cross-correlation is a (linear or nonlin- 
ear) function of the relative delay. The advanced signal 
will then be delayed by the amount computed from the 
cross-correlation. As mentioned before, the relative delay 
is not constant in time; therefore, to track the dynamics of 
the relative delay automatically, a closed loop is used. To 
close the loop, a number of components need to be added, 
e.g., a loop filter and an integrator. 

This article presents a closed-loop time-alignment sys- 
tem, which takes two baseband signals, estimates the rela- 
tive delay between the two signals using a correlator, and 
adjusts the time delay of the advanced signal through a 
loop filter and an integrator before combining the two 
baseband signals. The closed-loop time-alignment system, 
shown in Fig. 2, consists of a quadrature correlator, a loop 
filter, a unit delay, an integrator, a nearest integer device, 
and a decision-making device. 

Each component in the loop has its function, as briefly 
described below and in detail in the subsections. The 
cross-correlator determines the relative delay between two 
signals. The loop filter tracks the dynamics of the rel- 
ative delay and maintains the loop stability. The unit 
delay ( 2 -1 ) is required to close the loop. The integrator 
( z/(z — 1) transfer characteristic) combines the previously 
and currently computed relative delays. The nearest inte- 
ger device is needed because the delay can be adjusted only 
by an integer number of samples. Finally, the decision- 
making device decides which signal needs to be delayed. 


A similar time-alignment system called the real-time 
telemetry combiner (RTC) has been developed and ana- 
lyzed [2,4,5]; however, this closed-loop time-alignment sys- 
tem overcomes the following RTC limitations: 

(1) Since the length of the correlator in the RTC is fixed, 
the correlation feedback is not useful for low SNR 
and the loop has to be open. 

(2) Because cross-correlation uses the product 
of the sign bits of the two signals, it causes the 
quantization-noise power to increase from —58.8 to 
— 16.8 dB when compared with 8-bit quantization. 

(3) Since the delay can be adjusted only in one branch, 
the to-be-adjusted branch must be delayed relative 
to the other branch. However, the signals' ad- 
vance/delay relationship may change in time as the 
spacecraft moves and the Earth rotates. 

(4) The RTC has an analog feedback to adjust delays, 
which is not suitable for software or completely dig- 
ital hardware development. 

The closed-loop time-alignment system presented here, 
on the other hand, has the following features: 

(1) An adjustable correlation length to accommodate 
different SNRs. 

(2) Cross-correlation, which is a product of the 
floating-point samples of the signals, reduces the 
quantization-noise power by 42 dB with 8-bit. quan- 
tization. 

(3) A decision-making device that automatically delays 
the advanced signal. 

(4) Discrete-delay adjustments suitable for software or 
completely digital hardware development. 

Although this time-alignment system can adjust delays 
only to the precision of an integer number of samples, ac- 
curacy is not a problem when the subcarrier is a perfect 
square wave and the data transitions are instantaneous. 
In practice, accuracy is not a problem if the signal is not 
filtered or the sample rate is much higher than the subcar- 
rier frequency. In the case of a filtered signal sampled at 
a low rate, an added interpolation on the input waveforms 
can provide better resolution. 

Due to all these differences, the RTC analysis does not 
apply to the closed-loop time-alignment system proposed 
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here; therefore, an analysis of the new system is needed. 
Because of the nonlinear device that takes the nearest in- 
teger in the loop, linear system theory, or analysis in the 
frequency domain or Z-domain, does not apply. The re- 
maining choice is to use time-domain analysis for nonlin- 
ear systems even though it is not as straightforward as 
frequency-domain analysis for linear systems. 

A system analysis will be given later, following a more 
detailed description of each of the components in the 
closed-loop time-alignment system. Simulation results and 
comparisons with the theoretical calculations will also be 
presented. 

II. System Description 

A block diagram of the closed-loop time-alignment sys- 
tem is shown in Fig. 2, and an equivalent diagram of the 
closed loop is shown in Fig. 3. In this section, each of 
the components and their functions are described in de- 
tail. The components include a cross-correlator, loop fil- 
ter, unit delay, integrator, a device that takes the nearest 
integer, and a device that decides which signal to delay. 

A. Cross-Correlator 

To find the proportionality between cross-correlator 
output and relative delay, consider the noise-free case. 
The cross-correlator, shown in Fig. 4, computes the rela- 
tive delay by taking the difference between the correlation 
of Signal 1, si(k), with Signal 2 delayed by Nd samples, 
$ 2 (& — Ad), and the correlation of Signal 2, S 2 (k), with 
Signal 1 delayed by Nd samples, s\(k — Nd)] the result is 
divided by two. The correlator output becomes 

c(m) = t [ Sl (k) S2 (k - N d ) - »i(k - N d )s 2 (k )] (3) 

where s\(k) and So(k) are defined in Eq. (2). 

This cross-correlator is similar to that in RTC [2,4,5], 
except the cross-correlation is floating point instead of sign 
bit, and discrete instead of continuous. By using correla- 
tions of 8-bit quantized samples, the quantization-noise 
power reduces from —16.8 dB of the single-bit correlation 
to -58.8 dB. 

The analysis of this cross-correlator is similar to that in 
[2] except the cross-correlator is discrete instead of contin- 
uous. The proportionality of the cross-correlator output 
and the relative delay is given in Appendix A. Assuming 
that the signal power is normalized to 1, the correlator 
output in the linear region (|m| < Nd) is 


_ 4/jc Rsym _ 
c(m ) = m 

Is 

— am (4) 

where f sc is the subcarrier frequency, R S ym is the symbol 
rate, and f s is the sampling rate. Dividing the correlator 
output by a, one obtains the relative delay m. 

Appendix A shows that when the artificial delay, Nd, 
is one-quarter of the number of samples in a subcarrier 
period, N sc , the linear region is the longest, which means 
that larger delays can be tracked. Hence, the correlator is 
called a quadrature correlator. 

In a noisy environment, this output is further averaged 
over N correlator-output samples to reduce the noise ef- 
fect. 

B. Loop Filter 

To close the loop, one needs to add a loop filter for 
dynamic tracking and loop stability. A higher order loop 
filter does better dynamic tracking. In this case, how- 
ever, the relative delay is considered to have low dynamics; 
hence, only a first-order loop filter is used. That is, the 
filter transfer function is 

F(z) = K 

The choice of K influences the loop stability, loop SNR, 
and loop bandwidth. These issues will be discussed further 
in Section III. 

C. Integrator, Unit Delay, and Decision Making 

The integrator combines past and current delay compu- 
tations; the unit delay is required to close the loop; and the 
decision-making device is basically a two-way switch that 
switches one way when input is positive and switches the 
other way when input is negative. This decision-making 
device allows the system to delay whichever signal is ad- 
vanced. 

III. System Analysis 

To characterize the system quantitatively, a number 
of parameters need to be determined. Among them are 
the filter parameters for loop stability, the loop SNR, the 
loop bandwidth, and the SNR degradation due to mis- 
alignment. Except for loop stability, which depends on 
the filter parameters only, other parameters such as loop 
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SNR, loop bandwidth, and SNR degradation depend on 
the loop-output noise variance. This noise variance is an- 
alyzed following a brief analysis of loop stability. 

A. Loop Stability 

For stability analysis, the noise input in Fig. 3 is ig- 
nored because loop stability depends solely on the loop 
parameters and is input independent. Assume that the 
relative true delay is an integer multiple of the sample 
period. This assumption makes the nonlinear nearest in- 
teger device negligible in a noise-free environment. Note 
that when the delay is not an integer multiple of a sample 
period, the nearest integer device can be modeled as an 
additive quantization noise source of finite magnitude and 
thus can be ignored for the stability analysis. 

Consider a first-order loop filter, F(z) — K . The open- 
loop transfer function is 


and the closed-loop transfer function is 


1. Noise Variance at Correlator Output. To de- 
termine the noise variance at the output of the correla- 
tor shown in Fig. 4, assume that the samples of the two 
inputs to the correlator consist of the sum of baseband 
signals, s\(k) and s 2 (fc) as defined in Eq. (2), and inde- 
pendent samples of additive white Gaussian noise, n\{k) 
and n' 2 (fc), both with zero mean and variance <rf n . Denote 
the Arth sample of the inputs of the quadrature correlator 
as 

r .'( fc ) = «.•(*) + n'i{k) ( 6 ) 

where i — 1,2. 

The output of the correlator before “sum and dump” is 

«(*) = \V[(k)r' 2 (k - N d ) - r[(k - N d )r' 2 (k)} (7) 

The sum-and-dump filter takes the average of its input 
over N samples, and its ouput is then divided by the slope 
a to produce an output sample 


H{z) 


Gjz) 

1 +G(z) 




1 

aN 


( m + 1 )N 

E 


k — mN + 1 


q(k) 


( 8 ) 


A root-locus plot of this system with a first-order loop 
filter is shown in Fig. 5. From this plot, it can be observed 
that when K > 1, the root, locus is outside the unit circle; 
thus, the loop is unstable. For 0.25 < K < 1, the system 
is underdamped; for K = 0.25, it is critically damped; 
and for 0 < K < 0.25, it is overdamped. Besides affect- 
ing system stability, the choice of K also affects the noise 
performance in the delay estimates, as will be discussed. 

B. Noise Analysis 

To analyze the noise, we first determine the amount of 
noise (mean and variance) at the output of the quadrature 
correlator in Fig. 4. This amount of noise is added at the 
output of the correlator, which acts as a unit delay in 
the closed loop, as shown in Fig. 6. The noise statistics 
at the correlator output are given next, followed by the 
noise statistics of the closed loop. Because of the nonlinear 
device of taking the nearest integer, the Z-domain linear 
system theory cannot be applied here, and time-domain 
analysis is used instead. 


where a is defined in Eq. (4). 

The problem here is to find the variance of w m . With 
the assumption that {u; m } is a stationary sequence and 
the noise samples at the inputs of the correlator are in- 
dependent and identically distributed with zero mean, the 
variance of w m is found to be 

~ (2a ) 2 N 4 ^n ^2 + 4<r“ n2 Pi] (9) 

which is independent of m, and where P iy i = 1,2 is the 
2 th signal power, and of ni is the 2 th noise variance at each 
correlator input. See Appendix B. 

2. Noise Variance at Loop Output. An equivalent 
diagram of the closed loop is shown in Fig. 6, where noise 
is the input and the output corresponds to the location 
of the computed delay. Denote the k th noise samples at 
the input of the loop as in*., those at the output of the 
integrator as xjt, and those after taking the nearest integer 
as y*. The goal here is to find the variance at the output 
after the nearest integer device, that is, the variance of y. 
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The difference equation that describes the closed loop 
in Fig. 6 is 


where a 2 T is the linear loop-output noise variance in units 
of squared radians. 


Xk - Xk-i — K(wt-i - Vk- 2 ) ( 10 ) 

Assuming that w and x are stationary Gaussian random 
processes with zero mean, then y also has zero mean, and 
the variances of and y, denoted by <r~ , <r~, and <T“ , 

are found to have the following relationship (Appendix C): 


However, the loop-output noise variance, has units 
of squared number of samples. Converting the units to 
squared radians by equating a symbol duration to a cycle 
or 27 t radians, one has 


.> / 27 rR. 


tym \ 2 

V 


(16) 


2 A o .,•> o t-°2 

- h = h ' <T w 


(11) where f s is the sampling rate, and R 3ym is the data rate. 
The loop SNR can be expressed as 


The relationship between the variances of x and y is 

ff 2 = £p t * 2 (12) 


where 


Pk = 


-r 

j 


1 


x/2? 


:e -* 2 /(2 ° D dx 


(13) 


fs 


27T R S ym J 


(17) 


The one-sided noise equivalent bandwidth for a linear 
system is defined as [4] 


Bl 2 T l 72 


(18) 


One can solve Eqs. (11), (12), and (13), for <t~ for a given 
input noise variance cr^. Note that in the absence of the 
nearest integer device, = cr“, and Eq. (11) becomes 


where 


h = 


= Grr— I H(z)H(z~ l )-dz 

\H( Z )\ 2 27fj J\ z \ = i ’z 


2-K 2 -I< 2 
A' 2 + K ** 


*> 



(14) 


which is exactly the same as Eq. (D-3) in Appendix D, in 
the linear system case. 

Also note that for a large input noise variance (<r^ > 
0.5, observed from simulations), the noise generated by 
the nonlinear nearest integer device is approximated by 
a quantization noise with a uniform distribution and a 
variance of one-twelfth, or 


n 




H{z) is the closed-loop transfer function, Tl = NT S is 
the correlation time, N is the number of samples in the 
correlation, and T s is the original sample period. The loop 
input and output variances are related through / 2 , that is, 
var out = / 2 ^arj n . Note that uar m and var out are not the 
same as and al nt used previously. 

In a nonlinear system, however, Z-domain analysis is no 
longer valid; thus, there is not an expression for H(z), nor 
for / 2 . Therefore, the loop bandwidth cannot be defined 
the same way as for linear systems. However, an expression 
for the relation between the input and output variances is 
available in Eq. (11); hence, the one-sided noise equivalent 
bandwidth is defined in this case as 


C. Loop SNR and Loop Bandwidth 

The loop SNR is usually defined as [1] 


Bl = 


1 var out 
2 Tl var in 



(15) 




(19) 
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where a * and cr^ have been defined in Sections III.B.l and 
III.B.2. 

D. SNR Degradation Due to the Loop 

The magnitude reduction of the signal at the output of 
the sum-and-dump filter due to time misalignment in the 
baseband-combining scheme is (see Appendix E) 



where r, is the relative delay over a symbol period be- 
tween the ith antenna and the reference antenna. Note 
that the expression in Eq. (20), derived in Appendix E, 
differs from the stated expression in Eq. 39 in [1] for the 
baseband-combining assembly containing the RTC. The 
first and second moments of the reduction function of the 
time-alignment loop are 



where i — 1,2. Note that C Tl = 1, and = 1, which 
implies a T ~ 0. 

Assuming perfect carrier and subcarrier tracking and 
symbol synchronization, the degradation due solely to the 
closed-loop time-alignment system is 

D ta ~ 10 log 10 

E 7i c ?, + E 7i7j ('t, C t , 

i = 1 '■}<'*} /oox 


where 7 , = ( Pi/N 0 ,)/(Pi/N 0l ), and T = Ef=i7t. with 
N a being the number of antennas in the array, which is 2 
in this case. 


Having analyzed and defined all the necessary param- 
eters, we now can check theoretical versus simulation re- 
sults. 

IV. Simulation Results 

To verify the theoretical results, the P t /N 0 at the in- 
put of the cross correlator is fixed at 45 dB-Hz. The loop 
gain K is set at 0.25, and the sample rate at 500 kHz. 
The subcarrier frequency is 5 kHz, and the symbol rate 
is 1000 symbols/sec. There is no carrier in the simulation 
since the carrier is tracked and removed before baseband 
combining. The cross-correlation length N (or loop band- 
width Bi) is varied, and the loop SNR, p, is obtained from 
the theory and the simulations. The results are shown 
in Table 1 and the comparison of the results is shown in 
Fig. 7. 

The simulation results also show that when the noise 
standard deviation at the input of the closed loop, a w , is 
larger than the one-sided length of the linear region, Aj, 
{&W > Nd) the loop fails to track. This can be observed in 
Fig. 7 for p smaller than 17 dB. 

The SNR degradation due to the closed-loop time- 
alignment system is also measured through simulation and 
the results are very close to those obtained from Eq. (23), 
as shown in Table 2. The parameters in the simulations 
are as follows: Assume two identical antennas, each with 
Pt/N 0 ~ 23.54 dB-Hz. The modulation index is 70 deg. 
The sampling rate is 100,000 Hz; the carrier frequency is 
3 Hz; the subcarrier frequency is 1000 Hz; and the symbol 
rate is 200 symbols/sec. The time-alignment loop gain, 
A”, is 0.0625. Simulations are done for two cases, with 
a perfect time alignment and with the closed-loop time- 
alignment system. The difference in SNR degradation is 
due to the closed-loop time- alignment system. 

V. Conclusions 

The article gives a full description of the design of a 
closed-loop time-alignment system and a detailed descrip- 
tion of each of the elements in the system. It also provides 
an analysis of loop stability, noise effect, loop bandwidth, 
loop SNR in terms of the length of the correlation, and 
SNR degradation due to the loop. The closed-loop time- 
alignment system described here differs from the RTC in 
that this system is fully digital and is suitable for software 
and completely digital hardware implementation. The cor- 
relation is performed with floating-point samples of the 
to-be-combined signals rather than their sign bits. This 
causes the quantization-noise power to drop from —16.8 to 
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-58.8 dB in the correlation, assuming 8-bit quantization. 
The delay can be adjusted in either branch rather than 
a single branch only. The length of the cross-correlation 
can be adjusted for different levels of SNR. Because delays 
are adjusted digitally, an operation to take the nearest in- 
teger is needed, which results in a nonlinear system. An 
analysis of the nonlinear system is presented in the article. 
The noise analysis of this loop is significantly different from 
that of the RTC, due to the floating-point cross-correlation 
and the nonlinearity. The article also gives a well-matched 
comparison of theoretical and simulated results. The lin- 
ear region of the relationship between the cross-correlation 
and the delay is determined. When the standard deviation 
of the delay estimates exceeds a quarter of the subcarrier 
period, the relationship between the cross-correlation and 
the delay is no longer linear, which corresponds to about 
a 17-dB loop SNR. 

Although the time-alignment system described in this 
article aligns two signals, it can be easily expanded into 


the alignment of N a signals, where N a > 2, by aligning 
the signals pairwise in N layers, where N = ("log 2 No ] . It 
can be shown that the total number of alignments is A/a — 1, 
which is the same as if we chose the signal with the longest 
delay as the reference, and adjusted all the delays of the 
JV a — 1 signals one by one as in [1]. Two advantages of 
aligning the signals pairwise using this scheme are 

(1) No reference signal is needed as this reference may 
need to be changed in time. 

(2) The SNRs get stronger after each layer, thus en- 
abling the computation time to be reduced. 

The combining degradation still needs to be analyzed. 

Some simulation results show that this time-alignment 
system can be improved by low-pass filtering the to-be- 
correlated signals before the quadrature correlator. How- 
ever, further analysis is needed to completely characterize 
the system and the associated improvement. 
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Table 1. Performance of the time-alignment closed loop. 


N 

B l , Hz 

*l,(theo) 

<Ty(theo) 

Ptheoi dB 

Psimu } dB 

Psimu ~ Ptheoi dB 

25 

1821.6 

1303.6 

237.47 

14.26 

3.17 

-11.08 

40 

1157.1 

814.8 

150.85 

16.23 

9.82 

-6.41 

50 

926.6 

651.8 

120.79 

17.20 

17.13 

-0.07 

100 

463.7 

325.9 

60.45 

20.20 

20.13 

-0.07 

200 

232.2 

163.0 

30.28 

23.20 

23.32 

0.11 

400 

116.5 

81.5 

15.19 

26.20 

26.40 

0.20 

800 

58.7 

40.7 

7.64 

29.18 

29.13 

-0.05 

1600 

29.6 

20.4 

3.87 

32.14 

32.32 

0.18 


Table 2. SNR degradation due to the time-alignment closed loop. 


N 

B Li Hz 

Pi dB 

(theo) 

(Ty(simu) 

Dtheoi dB 

B simu > dB 

10,000 

0.172 

19.21 

75.97 

83.11 

-0.92 

-1.11 

50,000 

0.035 

26.18 

15.26 

15.48 

-0.39 

-0.38 

500,000 

0.0036 

35.96 

1.60 

2.45 

-0.14 

-0.13 

5,000,000 

4.58 X 10 -4 

44.96 

0.20 

0.00 

0.00 

0.00 
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Fig. 1. A diagram of baseband combining of two signals. 
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Fig. 2. A closed-loop time-alignment system. 
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Fig. 3. An equivalent diagram of the closed-loop time-alignment system. 
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Fig. 4. A diagram of the quadrature correlator. 



Fig. 5. Root locus for the first-order loop filter. 



















Appendix A 

Correlator Output Versus Relative Delay 


The following gives the proportionality between the 
cross-correlator output and the relative delay between the 
signals received from two antennas. Let N sc be the num- 
ber of samples in a subcarrier period, N sy the number of 
samples in a symbol period, and N d the number of sam- 
ples in the artificial delay. As stated before, a baseband 
signal consists of data on a square-wave subcarrier. As- 
suming that the data and the subcarrier are independent, 
the autocorrelation of the product of the data times the 
subcarrier equals the product of the autocorrelation of the 
data and that of the subcarrier. Considering the linear 
region |mj < N sc f 2, the data correlation function is 


Cdata{m) 


m 


N, 


sy 


(A-l) 


(m + N d ) 
N sy N sc / 4 


Finally, taking the difference of these two correlations and 
dividing the result by 2, one obtains an equation of the 
cross-correlator output in the linear region: 


c(m) 


4fsc - R 


sym 


fs 


m 


— am 


(A-6) 


where f s is the sampling frequency, f sc is the subcarrier 
frequency, and R sym is the symbol rate. 


and the subcarrier correlation function is 


c sc (m) = 1 - 


N sc /4 


It is clear that the longer the linear region, the larger 
the relative delay that can be tracked. The linear re- 
gion is longest when the regions [~N sc /2 + N dl N d ] and 
(A-2) [-N d} N sc f 2 - N d } coincide. That is, when 


The correlation of data times subcarrier is 


N d — N $c / 4 


C(data)(sc) ( m ) 


|m| |m| ( rri 2 

AW4 ~~ W + N, y N. c /A 


(A-3) 


Then correlation curves are delayed or advanced by N d 
samples, which is equivalent to the correlations of 
- N d )s 2 (k) or si(k)s 2 (k - N d ). That is, in the region of 
[— N sc /2 + N d , TVd], the correlation function is 


C ( S U'lay^) 


1 , N d rn- N d 
AW 4 + N sy 


(m - Ntf_ 
N sy N sc /4 


(A-4) 


and in the region of [-N d , N, c /2 - N d \ 


m + N d m + N d 
W/4 N^y 


Hence, the cross-correlator is also called a quadrature cor- 
relator. 

Figure A-l shows a special case where the number of 
samples per symbol period is chosen to be N sy = 10, the 
number of samples per subcarrier period is N, c = 8, and 
the number of samples in the delay is N d = N,J 4 = 2 to 
demonstrate the correlation of random data, a square-wave 
subcarrier, and their product. Note that N sy and N sc are 
arbitrarily chosen with a condition that N sy > N, c . 

The first part of Fig. A-2 shows the correlation of s, (k — 
Nd.) and s 2 (k) = si(k-m), which is the correlation of Si(F) 
and s 2 (k) shifted to the right by N d samples. Similarly, 
the correlation ofs^F) and s 2 (k - N d ) = s d k -N d -m) is 
the correlation of «!(*) and s 2 (k) shifted to the left by N d 
samples. The second part of Fig. A-2 takes the difference 
between these two correlations and divides the result by 2, 
so that a linear region can be observed from —N d to N d , 
with N d being 2 in this case. Figure A-3 confirms that the 
maximum length of the linear region occurs when N d T 3 is 
chosen to be one-quarter of the subcarrier period. 
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Appendix B 

Noise Variance at the Correlator Output 


To find the noise variance at the correlator output, we 
initially find the first and second moments of the corre- 
lation. Assuming that the correlation {u; m } is stationary 
(and noting that the index m can be dropped), then its 
first moment is 

£{»> =4e% ( h) 

k-l 

= E{r[(ky 2 (k - N d ) - r[(k - 1 W d )r' 2 (k)} 

a V k - 1 


Denote the variances of the noises at the correlator inputs 
as of ni and <r? and the signal powers as Pi and P 2 , re- 
spectively. Expanding the product q(k)q(m), one obtains 


N N 


E{w2] = [2dW» a W 

' Jfc=l rn — \ 


-f 2<7i n2 P\f>m,k + 2at ni P 2 6 m]k 


~ a in 3 E { S l( k ) s l( k ~ 2Nd)}&m,k-N 4 


V k-l 


- aj n2 E{si(k + N d )si(k - N d )}6 m:k+Nd 


— si(Ar — N d )s 2 [k)} (B-l) ~ <T in 1 E {s 2 (k)s 2 (k 2N d )}6 mk -jv d 


The square of the first moment is 
{E{w})2 

' Jt = 1 m~l 

- Sl (k - N d )s 2 (k)} 

x E{s\(m)s 2 (m - N d ) 

- s i(m - N d )s 2 {m)} 

+ E{si(k)si(m)s 2 {k - iV<j)s 2 (m - N d )} 

- E{si(k)si(m - N d )s 2 {k - N d )s 2 (m)} 

- E{si(k — N d )si(m)s 2 (k)s 2 {m - N d )} 

+ E{s,(k - N d )si(m - N d )s 2 (k) 

x s 2 (m)} (B-2) 

The second moment of w is defined as 

E{w2} = S ^ E{q(k)q(m)} (B-3) 

^ k = l m = 1 


- a f ni E{s 2 (k + N d )s 2 (k - N d )}6 mtk +N d 
+ E{si(k)si(m)s 2 (k - 7'/ d )s 2 (m - N d )} 

- E{si(k)si(m - N d )s 2 (k - iVd)« 2 (in)} 

- ^{si (X: - N d )s 1 (m)s2(k)s2(m - N d )} 

+ jE7{s! (X r - N d )si(m - N d ) 

x s 2 (X)s 2 (m)}] (B-4) 

For a quadrature correlator, the delay Nd is one-quarter 
of the subcarrier period; thus, two points of the same sig- 
nal that are 2 Nd samples apart are one-half of a period 
apart, which means that they always have the same mag- 
nitude and the opposite sign. Therefore, their product 
is —Pi, i = 1,2. Also note that the last four terms in 
Eq. (B-4) are the mean squared of w , which is subtracted 
from E{w 2 } to obtain the variance of w ; hence, the vari- 
ance of w is 

<?l = ( 2 0 )»Af ^" 1 g,,S " a + 4 ^ n3 Pl + 4a ' n ' (B ‘ 5) 
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Appendix C 

Noise Variance at the Loop Output 


The difference equation that describes the closed loop 
in Fig. 6 is 

Xk - Xk-i = K(w k -1 - yk-i) (C-l) 

Assuming that w and x are stationary zero-mean Gaus- 
sian processes, taking the expected value of Eq. (C-l), one 
obtains 

£{?/} = 0 

Then the second moments of the processes w, x ) and y are 
their variances denoted by cr^, cr “ , and <r^, respectively. 
Now taking the expected value on both sides of Eq. (C-l) 
squared, 

E{xl}-2E{ Xk x k - 1 }+E{xl_ l } = 

K\E{wl_ x ) + E{yl_ 2 }) (C-2) 

That is 

2a 2 x = K 2 al + K 2 a 2 y +2E{x k Xk-i} (C-3) 

The last term in Eq. (C-3) can be expressed in terms of 

the variance of x when multiplying Eq. (C-l) by x k -\ and 
taking the expected value. That is, 

E{x k x k - 1 } = E{xl_ x } + KE{w k -ix k -i} 

- I<E{xk-iyk - 2 } (C-4) 

To see that E{wk-iXk ~\ } = 0, one can multiply Eq. (C-l) 
by w k ■■ 


E{w n x „} = E{w n x n _ 1 } + KE{w n w n .i} 

- KE{w n y n _ 2 } = 0 (C-5) 

The last step in the above equation follows since the past 
outputs are independent of the future inputs when the sys- 
tem is causal and the noises at different times are indepen- 
dent, Going back to Eq. (C-4), since y k ~2 is the nearest 
integer of x k - 2 , it can be expressed through an c k -2 as 


Vk - 2 — x k - 2 + Q— 2 

where ejc -2 represents a quantization error. Substituting 
the above into Eq. (C-4), one has 

E{x k x k - 1 } ~ al- KE{x k -\x k - 2 ) 

- KE{x k -i(k~2} (C-6) 

The last term in Eq. (C-6) is zero since the past output 
Cjt _2 is independent of the future input x k -i- Applying 

the stationary property of x and rearranging Eq. (C-6), 
one obtains 

E{x k x k -i} = Y^7K a * (C-7) 

Substituting Eq. (C-7) into Eq. (C-3), one has 

- K^l = A VJ (08) 
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Appendix D 

Noise Variance at Loop Output for Linear Cases 


Without taking the nearest integer, the linear system 
theory should be applied. The variance at the output of 
the loop y equals the variance of the delay estimates due 
to the noise n, and it is 

= |;/( i )| 2< ^ /2 


where 




Using the formula given in [3], h is 

7 K 2 + K 
2 ” 2 - K - K 2 


The variance of y is 


Using the first-order loop filter, the closed-loop transfer 
function is 


A 2 T A o 
2 - A' - A' 2 


(D-l) 


(D-2) 


(D-3) 
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Appendix E 

Signal Amplitude Reduction Due to Misalignment 


Let N sc , N, y , and JV rci be the number of samples per q _ 

subcarrier period, symbol duration, and in the relative de- A ideal 

lay. First, consider the ideal case where the signals are 

perfectly aligned, and the combined signal amplitude at 2^ sy ^ rd ^ ^ rd 

the output of the sum-and-dump filter is “ jV 3C N sy 2 N sy 


Aideal — 27V$y 



Then consider adding two misaligned signals. The signal 
amplitude at the output of the sum-and-dump filter be- 
comes 



(E-2) 



1 N sc 
+ - 2 — 

2 2 


where f sc is the subcarrier frequency, R sym is the symbol 
rate, and 


+ 




(E-l) 


Ti 


Nrd 

Nsy 


The amplitude reduction of the misaligned signal at the is the relative delay between the ith signal and the refer- 
output of the sum-and-dump filter is ence signal in terms of a symbol period. 
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The real-time array feed compensation system is currently being evaluated at 
DSS 13. This system recovers signal-to-noise ratio (SNR) loss due to mechanical 
antenna deformations by using an array of seven Ka-band (33.7-GHz) horns to 
collect the defocused signal fields . The received signals are downconverted and 
digitized , in-phase and quadrature samples are generated, and combining weights 
are applied before the samples are recombined. It is shown that when optimum 
combining weights are employed , the SNR of the combined signal approaches the 
sum of the channel SNRs. The optimum combining weights are estimated directly 
from the signals in each channel by the Real-Time Block II (RTB2) correlator; since 
it was designed for very-long-baseline interferometer ( VLB! ) applications, it can 
process broadband signals as well as tones to extract the required weight estimates, 
The estimation algorithms for the optimum combining weights are described for 
tones and broadband sources. Data recorded in correlator output files can also 
be used off-line to estimate combiner performance by estimating the SNR in each 
channel , which was done for data taken during a Jupiter track at DSS 13. 


I. Introduction 

The advantages of array feed combining for recovering 
signal-to-noise ratio (SNR) lost to mechanical deforma- 
tions of large receiving antennas have been described in 
previous articles [1,2]. Typically, SNR losses become sig- 
nificant when carrier wavelengths smaller than the design 
tolerance of the reflector are employed. This is the case in 
the DSN, where Ka-band (33.7-GHz) reception is contem- 


plated with the large 34- and 70-m antennas, whose sur- 
faces are subject to considerable deformation from gravity 
and wind. 

The array feed compensation system currently being 
evaluated at DSS 13 has been designed to recover SNR 
losses at both low and high elevations, where losses are 
most severe. The idea is to collect some of the deflected 
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